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A simple model for the angular dependent interaction be- 
tween Ceo molecules in face centered cubic lattice is proposed 
and analyzed by use of the rigorous bifurcation approach. The 
quantitative results for the orientational phase transition and 
the characteristics of the ordered phase are in good agreement 
with the experimental data. 



The orientational phase transition in solid 6*60 is of 
much current interest. The orientational ordering has 
been a subject of extensive experimental investigations 
(see, e.g., jj]-5| ); some theoretical researches were per- 
formed, too 1 7 [KJ. However, ah initio calculations fail 
to reproduce the experimental results. 

In this Letter we develop a simple model for the an- 
gular dependence of the intermolccular potential in solid 
Ceo- The model is based on the ideas of preferred ori- 
entations due to David et al. 0] and to Lapinskas et 
al. [ pT[ , and on the maximal exploit of symmetry con- 
siderations. We apply to this model interaction the rig- 
orous approach based on the Lyapunov-Schmidt theory 
of bifurcation of solutions of nonlinear integral equations 
and obtain quantitative results for the transition tem- 
perature and the distribution of molecular orientations 
in the ordered phase. These results occur to reproduce 
the experimental data. 

As is established in a number of experiments Cqq crys- 
tallizes in a face centered cubic (foe) structure. At ambi- 
ent temperature the molecules rotate almost freely with 
centers on the fee lattice sites, so that the space group is 
Fm3m (see, e.g., g]). When the temperature decreases 
to Ts ~ 260A" the first order orientational phase transi- 
tion takes place: the sites of the initial fee latice become 
divided between four simple cubic sublattices (see fig.l) 



with own preferable molecular orientation in each sub- 
lattice. The broken symmetry space group is Pa3. 

Moreover, the neutron-diffraction experiments Q| 
have shown that the orientations in the ordered state 
are so that the electron-rich regions (the interpentagon 
double bonds) face the electron-deficient regions of the 
neighboring 6*60 molecule: the centers of pentagons or 
the centers of hexagons. It was shown jj],[5||j| that the 
ratio of the number of molecules in those two states is 
about 60 : 40 at the phase transition temperature and 
increases when the temperature decreases. This remain- 
ing orientational disorder is usually believed to cause the 
orientational glass transition at Tq « 90K now confirmed 
by various experimental technics (see, e.g., []l4|]). These 
two minima of the intermolccular angle dependent energy 
were obtained by numerical calculations and were shown 
to be much lower than the energies of other mutual ori- 
entations of the pair of molecules (see, e.g., P,P, pl||l5| ]). 
In those calculations the previously obtained charge dis- 
tribution for the isolated 6*60 molecule Jl(| was taken 
into account. Usually recent calculations use the inter- 
molecular potential of Sprik et al. fll7|| : a sum of 6-12 
and Coulomb interactions between 60 atoms C and 30 
double-bond centers D and between each other: 
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Here e = 1.293meV,acc = 3.4A, ocd = 3.5A, <7dd = 
3.6A, q D = -0.35e,<7c = -to/2. 

Rigorously speaking we are interested in the angular 
part of this complicated interaction represented in terms 
of multipole-multipole interaction of point-like multi- 
poles on the sites of rigid fee lattice with coefficients to 
be calculated from ([!]) . The general form of this angular 
part is 



C* (0Jij )ui v {ui)ui T (U!j ) , 



(2) 



Fig.l. 



with I = 6, 10, 12, 16, 18, ... due to the icosahedral molec- 
ular symmetry //,. In the Eq. (0) uii are the angles de- 
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scribing the orientation of the molecule on site i, for ex- 
ample, Euler angles and u\ v - some kind of harmonics. 
However, we simplify the problem and develop a model 
orientational interaction. As the angular dependent in- 
teraction is rather short-ranged we can restrict ourselves 
by the nearest-neighbor interactions. 

We follow the main ideas of the papers 
use the restricted number of allowed orientations instead 
of free continuous rotations. Let us take into account 
in the energy (|J) only the orientations with pentagons, 
hexagons or double bonds directed towards 12 nearest 
neighbors in fee lattice. The Cqo molecule is constructed 
in such a way that if 6 of its 12 pentagons (or 6 of its 
20 hexagons) face 6 nearest neighbors double bonds (P 
and H states of Lapinskas et al. ifTllO , then 6 of its 30 
interpentagon double bonds face the remaining 6 nearest 
neighbors. Now the energy matrix elements can take 
only three values: Jo - the energy of the general mutual 
position, Jp - pentagon versus double bond and Jh - 
hexagon versus double bond. These energies in our model 
can be compared with those calculated in [j§|,|9| |l3] , |l5| ] as 
functions of the angular displacements of the molecule at 
(0,0,0). Following Q, and putting J = we obtain 
from the fig.2(b) of the paper g J P = -300X and 
Jh = —110K. Now we leave the paper and follow 
our own way. 

The energy matrix elements Jp and Jh connect the 
states of molecules only in the allowed orientations. So, 
only allowed linear combinations of ui v enter the Eq. (|^). 
The theoretical curve in || makes no difference between 
the number I of harmonics and describes the effect of all 
of them. So, in the framework of our model calculation 
it is possible to build up the allowed functions using only 
the harmonics with 1 — 6: we need only their transforma- 
tion properties. We restrict ourselves to I = 6, however 
the coefficients Jp and Jh are not some of C® T given in 
(||) but effectively take into account higher order terms. 
Let us construct the functions Pi(uj) and Hi(u) explicitly 
in terms of cubic harmonics K m = K§. m , m — 1, 2, 13 
(see, e.g., Jh|). All functions Pi and Hi are the sums 
of K m , invariant under the icosahedral symmetry of the 
molecule (i.e. belonging to the A ig representation of the 
icosahedral group Ih) if icosahedrons are naturally ori- 
ented in one of 8 properly chosen coordinate systems. 
The states Pi(Hi) have 6 pentagons (hexagons) and 6 
double bonds directed towards 12 nearest neighbors along 
different [100] axes. Piioj) describes the molecule ro- 
tated from the standard orientation B (following 18|) 
about [111] axis through the angle 97.76125°. The angle 
for Hi\w) is 37.76125°. The functions P 2 {uj),P 3 (uj) and 
Pa(lo) (or H 2 , H 3 , H4) are obtained from Pi(u>) (Hi) by 
subsequent counter-clockwise rotations of the molecule 
by 90° around z axes. 

If written in the standard coordinate frame with Carte- 
sian axes along the cube sides these functions have the 
following explicit form: 

Pi(w) = a P K x {w) + (3 P [K 8 (u) + K 9 (u>) + K w {u)\ 



+ jp[Ku(w) + K 12 (u>) + K 13 (oj)], 
P 2 (lo) = apKtiu) + Pp[-K s (lu) + K 9 {u) - K 10 (u)] 

+ jp[-K n (u) + K 12 (uj) - Ki 3 (w)], 
P 3 H = ctpKiiu) + p P [Ka(w) ~ K 9 (lu) - K w (co)] 

+ j P [K n (u) - K 12 (uj) - Kisiu)], 
Pi(u) = a P Kx{u)) + Pp[~K s (lo) - K 9 (u) + K w (to)} 

+ 7p[-K u (u) - Ki 2 {w) + K 13 (u)], (3) 

with a P = -0.38866; (3 P = 0.31486; j P = -0.42877. The 
functions Hi(u>) have the same form as Pi(ui) but with the 
coefficients a H = 0.46588; j3 H = 0.37740; 7^ = 0.34432. 
The functions are normalized to unity. 

Let us now treat our model by use of bifurcation 
approach in the mean-field approximation. As is well 
known the mean-field approach often brings one to the 
formulation of the broken space symmetry problem in 
terms of the bifurcation of solutions of nonlinear integral 
equations for distribution functions (see, e.g., the review 
Jl9|). In particular, the bifurcation approach was used in 
the case of orientational phase transitions in molecular 
crystals in the Refs. pfl|-p5[] etc. The simplified version 
was originally developed by James and Keenan for solid 
methane p6| and by Michel, Copley and Neumann Q] 
and by Heid ||l 0f| for solid Cqq . We shall follow our papers 
on hydrogen |20|- p2] , p5[ . 

In the mean-field approximation from the first equa- 
tion of BBGKY hierarchy for the orientational distribu- 
tion functions or by minimizing the orientational free en- 
ergy one can obtain the following nonlinear integral equa- 
tion §i: 



(4) 



gi(u>i) — lnf^gr 1 ^], fi(u>i) - one-particle orientational dis- 
tribution function for a molecule on i-th lattice site, the 
constants Gi are the normalization constants. 

In our case of solid Cqq where there are four sublat- 
tices (see fig.l) and four kinds of unknown distribution 
functions we obtain from (|]) the following system of four 
nonlinear integral equations: 



+ A(uj,uj')G z e 93 ^' ) 

+ A(uj,Lj')G4e 9i( - UJ ' ) 

+ A(u,u')Gie 9l( - u ') 

+ A(u,w')G 2 e 3 ^ 
(5) 



5 i(w) + A / doj'[B(oj,oj')G 2 e 92 ^ 
+ D(u,u')G 4 e 94 ^">} =0, 

g 2 (uj) + \ J duj'[B(Lj,oj')Gxe 91 ^ 
+ D{uj,Lu')G 3 e 93 ^">] =0, 

g 3 {cj) + \ J dcj'[B(u,u')G 4 e 9i ^ 

+ D(u,uj')G 2 e 92 ^] =0, 
34 M + A J dui'[B(ui,ui')G 3 e 93 ^ 
+ D{uj,Lu')G ie 9l ^">} =0. 
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Here A = 1/T, A(u,u'), P(w,u/), D(u,lu') are the sums 
of interactions over nearest neighbors in the sublattices 
A, B and D (see fig.l), respectively. For example, the 
sum in the plain perpendicular to the x axis can be writ- 
ten explicitly in the form: 



D(lu,lu' 
2{[(J\(w 

x[P 2 (u;' 

+ [P 2 (lu 
x[(Pi(u/ 
+[(P 2 (w 

x[Pi(w' 
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x[(P 2 (w' 



P 4 (w)]Jp + (J?i(w) + Jff] 

■ P 3 (J) + H 2 {u') + H 3 (u')] 

■ p 3 h + fr 2 (w) + h 3 (u)] 

• P 4 (o/))J P + + fl4(w'))^] 

P 3 M) J P + (P 2 (u;) + ff 3 (w))4] 
•P 4 (a/)+Pi(a/)+P 4 (w')] 
P 4 (w)+#i(o;)+P 4 (u;)] 
P 3 (a/))J P + (H 2 (w') + fl3(w'))J H ]} (6) 

and analogously for two other sublattices. 

The equations (|^) are well known Hammerstein equa- 
tions ^7|. In the case of finite domain of integration 
when the fixed point principle is valid there exists de- 
tailed theory for such equations (see pi]). We use the 
standard methods (see, e.g. [^9|). At high temperature 
the system (j^) has only trivial solution ^(o;,) = 0, corre- 
sponding to the orientationally disordered phase. At the 
bifurcation points A Q new solutions with broken symme- 
try appear ( X a > 0). For A = A Q (1 + /i) the functions 
gf (tOi) can be written as series in integer or fractional 
powers of fj,. These powers are defined by the bifurcation 
equation (see [^9| ) corresponding to the system (|J) . In 
our case we have 

gi{u) = fJ,hi(u>) + ^ 2 Xi(u) + ... 

because among the integrals / dujK mi K m2 K m3 with m = 
8, 13 there are some which are not equal to zero. This 
means the first order phase transition [p0|-p2[ . 

The bifurcation points are the eigenvalues A Q of the 
linearized system corresponding to (J5|) : 

hi(u) + — I du'[B(u,cf)k2(t^+A(u,u/)h 3 (u r ) 
An J 

+ D(uj,uj')h 4 (u')} = 0. 



The functions hi can be written in the form 



(7) 



(8) 



have truncated the interaction and reduced the problem. 
Nevertheless, there remain still two quantitative charac- 
teristics we aim to obtain: the bifurcation temperature 
and the relation between the weights of P and H func- 
tions in the solution. 

Using (||) it is easy to rewrite the system (0) as the 
system of linear algebraic equations for the coefficients 
h% with i = 1,2,3,4 and v = 1,8,9, 13. 

Using the explicit form of the matrices A, P, D it is 
easy to obtain the only nonzero elements: 



A 1A = P M = A,! 
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One can write the elements u, v, z and w in terms of the 
coefficients ap, (3p,^jp,aH, (3h,1h and energies Jp,Jn 
and obtain the following values: u — 32 • 5.046, v = 32 • 
94.127, z = 32 • 7.665, w = -32 ■ 37.155. 

The determinant of the algebraic system is factorized 
in 2 x 2 determinants, so that the eigenvalues A Q can 
be easily obtained. Among the values X a there are two 
positive values. The first one Ai = An ju corresponds to 
the solution proportional to K\ and is of no interest now. 
The second is the positive solution of the equation 



1 A / \ A 2 

l -^ {v + z) + TAny 



(vz 



w 



0, 



(9) 



namely A 6 = 0.00364X- 1 or T b = 275K. The corre- 
sponding nontrivial eigenfunctions have h\ = 0, and the 
other coefficients h^ are subject to some constraints. If 
we add the condition for the functions hi(ui) to transform 
one into another under the action of the cubic group rota- 
tion elements which leave the fee lattice invariant, then 
only three of the coefficients remain to be independent 
and the functions hi can be written in the following form: 

= aPi(w) + bH t (u) + cKx{u)), 
h 2 {uj) = aP 3 {uj) + bH 3 (Lu) + cK x {w), 
h 3 (uj) = aP A (uj) + bH A {uj) + cKi(u), 
h 4 (uj) = aP 2 (uj) + bH 2 (cu) + cKi(w), (10) 



so that the eigenvalues X a define not only the bifurca- 
tion temperatures but the relations between nonzero co- 
efficients h\ (that is the symmetry of the new phase), 
too. 

In the case of the full interaction one can obtain all pos- 
sible broken symmetry phases compatible with the ini- 
tial symmetry and the condition of positive temperature 
value (see, e.g., the case of hydrogen p0|,p2p^] ) . Now we 



aap + ban + c = 0, 



a(3 P + b(3 H = Q(aj P + bj H ), 



4"7T 



(11) 
(12) 
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Using the numerical value for Q we obtain immedi- 
ately: 

Pp = -4t = °- 608 ; Ph = t - r = 0.392, (13) 
a + b a + o 

so that not only the transition temperature, but also the 
ratio of the number of molecules in P and H states occur 
to coincide with the experimental data jl]|],|5| pp ■ Ph = 
60 : 40. 

To obtain the remaining unknown coefficient we use 
the equations of the second order in \x {giioj) = phi{u>) + 
p 2 Xi(ui); A = A c (l + p)). The system has the form: L x 
x = R, where L is the linear 4x4 operator L x h = 0. 
The system has nontrivial solutions for Xi if the right 
hand side is orthogonal to the solutions for hi obtained 
before. All 16 equations 



dujRi(ui)hj(uj) = 



(14) 



are identical due to th e sy mmetry of coefficients. Solv- 
ing equations (|ll"|), ( [l2| ) and (llif), we obtain finally 
a = —25.7; b = —16.6; c = —2.27. The minus sign 
means that the solution goes in the direction of higher 
temperatures (p = —t,t = (T — T&)/T(,). The solution 
has the turning point T t which is some Kelvins higher 
than Tjj. The actual first order phase transition obtained 
from the free energy behavior takes place between these 
two points. The details of this calculation are to be pub- 
lished elsewhere. 

To conclude, we developed a simple model for angle 
dependent interaction for Cqq molecules in the fee cu- 
bic lattice. We used rigorous analytic approach based on 
the Lyapunov-Schmidt theory of bifurcation of solutions 
of nonlinear integral equations to treat this model. As 
result we obtained the first order phase transition, the 
bifurcation temperature Tf, = 275K, the Pai symme- 
try of the ordered phase and the ratio p of the number 
of molecules with pentagon facing neighbor double bond 
to the number of molecules with hexagon facing neigh- 
bor double bond (|l3|) near the phase transition in good 
agreement with the experimental data. 
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